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Mme LUEEOSS crtk 
Bodificaticn and impr a 
(OPTSYS) used in tre study, GSSiG5, 260° “explicang 
Geeaisal Systems Contr 

Mbes 9¢rtinal sy &m 
emcees y Gali [ khef. 1] in 1971 to suprort his res*aich 
mieccednhyewing VIOL eircraft control systems The latest 
program mcdifications were mads py Walker {Ref. 2} and Lin 
(meres ct Stantord University, to OPTSYS 4 and CFTSYS 5 


lespectively. These fprograg versions performed quite satis- 
factcrily in the batch envirecnment, but exhibited varying 
Gegrees ci user hostailicy due to data input format requirze- 


ments and incomplete frogram documentation. 

; meemcsadanal intéert of this work was to adaot Walker's 
Peeeta-a Version of CETSYS to run interactively under Vu/CcHs 
ch the IEM 30333; hcewever, the extensive mcdificaticns of 
CPTSYSK rew represent a high-level incterective applicaticns 
software system capakle of integrated simulation, analysis, 
synthesis ard design of broad classes of optimal systems 
Gomrrel ~robléens. With OPISYSX users may now evaluate 
varicts specialized crtinal systems control applicatiors, 
relieved of the Eurden of lengthy efathematical fpregran 
develcement. 

It is assumed thet the reader/fuser is familiar with the 
tasic ccrcerts of Certrol Theory and Optimal Systems Design. 
The prektiea descripticns and frogram development fcilow the 
Geeianci¢cg y and symbol/fnaming conventions Or EZYSOr 
[Reft. 4}. 

An e€xplanaticn cf the tasic system of egquations, the 
terms and symbolcgy used, and a program overview including 


*he general methcds cf soluticn are presented first. 
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II. THE OPTSYSX COMPUTER PROGRA} 


a <2 4900 ‘eee TP Ok GES = 


ae INTECDUCTION 


Seleisk 25 a double-precision FORTRAN progrten erricving 
Deitel CCUcerCi tneory aralysis téchnigues. Ajkbouch the 
[ere Meboceoragdinally written to synthesize contrciiers ior 


a z 

BacmeOoL almerare {[ kof. 5], i nas been Bxteneively 
aC ean Le Conzrclier, PE Tens on eee 6 U leseos 
synthesis as well as transfer function and modai analysis on 
coher tyres of large, multi-variable systems or equaticns. 
jiagemercaram modifications described in this work nrew aliow 
Eavid frureracal computer analysis in a completely interac- 


tive mode. 


toes tern /TODEL DESCIIPTION 


CEISYSK treats a linear stationary system model: 
x = [Fx + (GJu + [GAMMA ]w eS 
cutput equation 
y = (H)x + (Dju (Zee) 
measurenent equation 
=f ix + [Djw + ¥ (23) 
where 
vu = centrel vector (m X 1) 
y = wkite meastrement noise vector (fp X 1) 


¥ = white process noise vector (q X 1) 


Meee taece vectcr (n xX 1) 
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Wares LR Ae ee kee yea CS USTLLY. (sya weet is Of 
MPePecieeex) , 1G | 1S the control distrisution matrix; 
(ear) 25 the process noise distribution matrix; {HH} is the 
Reacsurement distribution natrix; and and {D] may represent a 
Poem coWa-c distribution matrix cf either the process roise 
yect 


n value, and a -Cevacilarcs 


E(w) = 0 (2.4) 


(Cj) = Eww": (2.5) 


Meee y ec Onc has Zero mean value and a COvarliance matrix 
(Rj, where: 


E(v) = Q (25) 
and 
[k] = Evy") (2.7) 


Meemevedratic perfeimance (Or cost) index for the Linear 


guadradic regulator is the expected value cf; 
VS PAG IS ie Geis Hey eae (2.8) 


in the statistical steady-state, where [A] represents an 


GUteut ccst matrix (a weighting on the output variables); 


TZ 





Meet tescace Weighting. is desired, {H]}] is fEttresented 
€ 


Meee escehtlty ratrix {I}. 


Cw. EFRBCGEAM OUTPUT 


1. CrEen-Loor Eicensysten Caicuiaticens 





P25. SA Cut cue nC lUuGcesS t= 


Wikec, Se) ct 


C 

iene gjs Seuwuy tie USeCr Ter That, + 
cns béing modeled, 
c 


a 
TOE RCCLCII Mae nwOne Or c 





GuUlatOs CrCl ien, 


full state variatle feedback is assumed as, 
(C] = (B }{¢° }(S] (eee 
and 
ded OND (2.10) 


Mececonut rol cain {Cj 1S a matrix GE optimal gains 
which wirimize the ccst index expressed in equation (2.8). 

For optimal reguiator synthesis preblems, fprogram 
cutput includes the closed-locp eigenvalues and #¢igenvec- 
fomesseene Control gain {C]; the closed-loop dynamics matrix 
PeeGC}; anhe the steacy-state gain matrix {5S], where {S] is 


the steady-state solttion to the algebraic Riccati equaticn: 


S@retes iS - S(C}LE JiGis + [H}{a}[#} = 0 (2.11) 





= (FX + (K](z - (BX) (ae) 


x 2 


™ 
xX 1S the stete estimates 
Weemetcmamilecrax of filter gains, 


and tre state covariaérces is described by: 
= T AAT ‘ 
E(xx ) = Eixx ) + (P- (22-159 


Meemerectesr 4 «€6©6CaeiN Matrix ({K] of equation (2.12) is 


Chteaned frem the relationshir: 
oie fet H NR} (2. 14) 


ahnere (Fj) is the steady-state solution to the algetraic 


Ficcati é€quation: 


; U ree 
(FIE + PLP] - PCEICR CHIE + (GOS) = 0 (2. 15) 
representing the errcr covariance of the estimate <. The 
centr¢cl ccvariance is the expected value described by: 
E(uu') = (C7x[C ] (2. 16) 


For the Kalran filter/optimal estimator syntkesis 
PuombemseeeGrisYSX olteut includes the eigenvalues and eigen- 
vectcrs cf the optimal estimatcr (Kalman filter); the filter 


Gaumm= | KG the error covariance matrix [Pj]; the covariance 
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- Sm ae ma Pa fa. . 
Meee es OL che SStlitate fx]; the state covariance vatrixu 7; xi} 
described in equatior (2.13), where 

(ex = Fal xx (2.017) 
emamer=] COntrOl Covariance mattix {U] describsd in squetion 


feces. et aonary.r¢ 


s i 
Secmpbe=schecad e252 roct-mean-Square values of the state and 
i | 


g 
ecntrel covariance ratrices{[Xl and (U] described in ¢qua- 
coc (2517) and (2.16) respectively. 


Ce. SCLUTICN ALGCRITE# 


Cne cf the fundamental techniques necessary to q 
synthesis cf optimal contrel systems is the stea a 
Semple a2cn CE =he algetraic Riccati eyuation. MASS 12S 'a wen 
Remvadal task duc to the iterative nature of ¢t u 

Tre steady-state soluticn by any gq € 
necessitates tine ircrement selection no gre than seme 
Peaerwen cf the shortest period of the cl 


enditure cn 


Ss 
ie 
r 
sed-locp system; 
eposing a Significant computer solution 73 Pp 
e 


t- 
m2 8B O 
(D 
wD 


the usér as well as thé requiring an agcunt of 


© 
computer Storage carekility due to th expansion 
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Peeeecr invelved. PeEenoer 5 ssity cf a 


ct 
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(D 
a @ 
QO 
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WD 
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time-varying solution of these equations fer ptinal cren 
Ugeopmereercrc! Ct ¢€Stifation tregquires the inversion cf ar 
n X¥ n matrix for each time increment where an unsteady solu- 
wen 12 desired. 

Pepewe=aul and efftcient alteznete methcad cf scluticn 
was develcped by Bryscn and Hall [Ref. 7], cased cn eigen- 
vectcr deccmpositior of the ¢sigensystem of the constant 


coerficient EULER-LACKANGE equations. 
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Geechee OPtcamal regulater, these #quations tak< +h- 


toed: 

ee 7 7 ee ai a oes 

D4 = = Gs 3 ¢ 
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i = A -FT » 

Mee che Obco@al fiiter or eStinator, ths equaticnké ere 

Carers rcrn: 
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x = a 

® 

Tt O-l 
A] = BROW 


Gam*)*Gaa’ x 
am! 
IN 
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GCEIsysx, and all earlier program versicns, us¢ the 


rethed cf ¢igenvectcr decomposition or the fULER-LAGFANGE 
eguaticns descrired in ({Ref. 7] for guadratic sy = 
contrel systems. Prcegram caiculations are based cnr th2 Qk 


INSORwey Ret < faa). 


oO Ww 


aiectathe Ci Francis, modified by wil 
g 


MrecrEtamt fc€atuzes <i this tethcd of nsystem analysis 


2 
are its improved aACCULaCY, ard s2eSLhSenstcavity ce widely 
separéetec eigenvalues. 

A mcre detailed description of tne QR algorithm ard its 
numerical @pplicétiors to eigensystem analysis may Fe found 
adm { Ret. 6). 


Fe EROGFAM OVERVIEW 


Cres lo x and its 41 subreutines may obd¢ divided into 
Raines beste Categories; 
1) setup éerd computaticn sequencirg 
z) data infut 
ay calculaticn 
A krief and general description of the program accdules 
and eubrevtines suppcrting these basic categories ccrcludes 


tas) Sect icn. 


ie 





Gecacre. 


cata infgut is perrormed by ths 14 REAL__ 
Mhz Oba nes . Bo edeteaied deseraption  ~cr the operaticn of 
these etkreutines is also included in Chapter III. Internal 
BuecgetcrattoOn Of external data input is provided by 
motcleec SETUP. 


Bewpcealcubation Sequencing 





INNER furctions as a second MAIN perecram. 
iemwonmecer= tre data irfut/calculation sequences for the tyre 
cf preklier teing solved and performs numerous matrix calcu- 
eee DS . Ji is fren this subreutine that ail input and 


Calculaticrn sequences are order2ad and performed. 


Subroutines MINV, BALANC, ORTHES, ORTRAN, HCH, HQH2, 
Eaboan, CNCRH, and EREXIT perform the major calculation 
Pegmerecc= CL the "OR" algorithn. 


Bpmeebeccati Equation Calculations 


Subroutine FGAIN sefparate2s the eigenvalues and 
eigenvectcrs of the Euler-Lagrange equations by eigenvector 
de=cemececiticn. RGAIN and subroutine SCOV calculate the 
Eteady-scate soluticr cf tke Riccati equations fcr the 


Gememelwer Cr @stimatcr proklen. SUNG Omelne seGy COMEUTES 
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Simrcuene MGEE  Confutes the modal transtormeaticns 
required fcr modal aréelysis. 


Pewee Sobers FUNCT. On Calculations 


a a =D a ere aap «aR —_ Op ge 2. 2p a» mWeE-aee ae Se ee 


Qutines TIF, FOLES, ZEROS, RESID, ACOM:, CCOME 
BPowreeanc runctica scl 


i @} 


EOL teansSie=> LUnCTEac: GemeeuTaticss 
asscciated with Modal calculation sequences. Sablss 
ESDCAL ccmputes the [cwsr spectrai density of the cute 

ie 


ccntrcls cf a controlied systen. 


temeesace OUL CUL 


ebject time format. 


format screin viewinc cr a 


1& 





(Gea pa al) 


= Set PS SPS a eers wee -we? ae ae <=> ay ents can =n <n ane — Bt a a oe we i eee eo 


é ou 
€ us¢€r cperation. Experie 
Geave COMpUter ccamunicaticn cifers many advantas 
esearck and prekiem solving environment. Ihe G op Gite ul 
nity fer flexible and immeciate computer communicaticn, as 
weil &@S the ability to salect alternate soluticn methods, 
e Significant advartages to the user; advantages which are 
tnaveailartle in a batck processing envirenment. 

Althcugh previous varsions cf OPTSYS produced all the 
desired calculations inthe Fatcen environment, goo Cae lah jee 
format and data sequencing and naming conventicns were 
Gontusing tc many users. The user was burcaned with ths 
recessity cf verifying the correctness of input data fcrrat 
and prcgrat tlow control flag settings for éach program run, 
in crder tc ensure the desired caiculation sequence was 
properly performed. 

These requirements, ccmbined with incomplete prcgran 
Gocumenteticn, promceted a lack cf confidence in the results 
and discctraged many potential users srom é€ven attemptirg to 
us2 CETSYS. 


Ej. FERCGHAF LANGUAGE 


Mreevon 1S Erogremmed in FORTRAN iV, folicwing the 
converticns of IEM System/360/370 FORTRAN IV language. Very 
few pregram features tkave teen peen incorporated which are 
Peewee t=6 in ANSI Standard FCRTRAN. These subtle differ- 


Smeee aCtWwitn=tandirc, OPTISYSX has been compiled and run 





Opeer Ecth FORTRAN IVY (61) and FORTRAN 4 (ate Tec) 
nigem 


Meme llGTE Ch the is? 3033. eeu Sete Ghee et OO Lad 
emeegee os mee CR SS CL 2600 iires of taxt, it is corsicered 

Gompi= tcly forteflie ficm one cferating system to anctnec. 
Crect= pPreStumpticnh that mos= scientific and restarch 
perscnnel are familiar with FORTRAN language, program modi- 
:§ Mey be easily undertaken onca system oOperaticn is 


Cw. GENERAL FROGHAM ECDLIFICATICNS 


All cf the previcusly deéevelcped numerical calculatior 
séequerces cf OPTSYS were retained un-modified in OPTSYSX. 

Thosé program se€equences requiring the input cf diagonal 
ccst cr ccvariance matrix elements were deleted or modiried 
geen XhEECVe user flexibility in entering any d¢sired 
weighting matrix, diagonal cr ctherwiss. Ties mee lt. Caclon 


streamlined program cperaticn thrcugh 2limination cf several 


program fl¢ew control flags; reduced a measure of user urcer- 
tainty; anc further decreased the required degree of user 


Eeogiran fariliarity~-fromoting uninterrupted cperatica. 
ieee  “sunhOUutines LCepEéeséent variations om the prin- 
Gmercecr Sample, effective methods of interactive inptrt, 
coupled sith errer-ccrrecticn/recovery sequences. 
Sukrcutines aDCHAR, RDINT, 4&nd RDREAL were developed by 
woemetroaer and Cdr. F.D. Sullivan to accomodate varying 
Pome ret mat requirements; null-string entry protecticn was 
@ewelcped by Cdr. P.£. Sullivan. Thess program features are 


Geetmescg in greater detail later in this charter. 





Le INTERACTIVE FROGEAH DEVELCEMENT 


O 

FEograd input and calculaticr sequence. After careful anal- 
Moms  Ceey net, 2], &@ basic fEregran flo# c a 

Sere a jisred. Ptotvemenis £1ew Ccenrrol d 
Ezanching network was formulated whose path coul 
Mined through either kinary icgic or numerical se 

Gisecemepasie LEanching  catégqeries we 

from the various protiem dascription statements: 
1G) G] SKC E DS Secor er eNO! . 6} 
Pac egek=--—--= rf toutes © ora st 
moo linees—-—-—-={"Ineue tke Valle or... } 

From the viewrcinte ci free-format computer ccmruni- 
caticr, integer and real number input presented no Ssignifi- 
cant preotklers. FORTRAN ccmcriler language is written such 
that numerical data input (réal number or integer) is 
eepeetcces  thebeky legquiting cniy an INTEGER or REAL déta 
Meen statement Within the progtan. Once the data type has 
feen declared, the desired valués may then be input with a 
free-fcrmat READ (5,%*) statement. 

Cne note of caution ccencerning numerical data ingfut 
in tree-format deserves menticning: Although all FCETRAN 
compilers treat character string input és an iliegal data 
Sype CCUVErSion, mary Will automatically convert the inad- 
vertent character entry to a “zero" and continue executicn. 
EFrotecticn against inadvertent errors of this tyre is 
discussed léeter in this chafter. 

Logical input ("Yes" cr “No") poses a@ unigue precblen 
tc prcegriammers. tee Uucld LelertieOd Nomen CcelLpobat.oOn 15 06 
Bommere es tre uSer tc convert the legical answer intc an 
Demecersi.e., "Yes*® = 1, "No" = 2. These integers may then 


ke tread cirectly intc the program, determining prcgram flew. 


Z1 





Bletough this methed gay promote fpregranmeing ¢25s¢, 
2c EGuITEes tne ser CCunagoDpT So Wee ae echt, ta 
Oaee-oe- “Che Which itczrezses the possisility of ebnezzali 

ic 


é 
Meeecdccdyiptc OFISYS<. 
Peeslogicait Strings ("Yes end "No) are declared as 
Seeeece-r Strings in adata type statement with 
EHogte® Gr Subprcgrar. A fcrmat statement is alse included 
foeectec trecr2m Or subprogram utilizing the “A"-field tc 
Seecaty tle desired ckaracter field width. S REWER ES .22e— 
Bent 15 ther inccrporated in the specizic program or sukpzo- 
Saatesmm@edie-cly pricr to each logical string input peirt. 
This REWIND statement allows the input device (the términal 
Bae Case) tC read ths character strine in the 
em@iroz es £26e~-fcrmet deta input. Dhewenalbecrwer, I1¢re 
Meeker cr tras modification was established at Al, alicwing 
streamlined operaticr with the user typing either "¥" er 
mee rer an affirmative rerly; "N" or “No for a negative 


ECDL Y « 
Z. General input Sequencing Requirements 


Peemedete @o5ue tO OFTSYSX 1S ain matrix or vector 
EOE Da = « This data input must be correlated in acccrdance 
with the preblem description and then properly sequenced in 
Gate=erece etre prcogran tc pericrm the desired calculaticns. 

fememorrginal and medified OPTSYS programs [Ref. 2] 
Meneame a yeregusred not only fFroblem description knowledges 


Eu cemelete user faniliarity with the detailed calculaticn 


séquerce cf the frogran. The latter point was considered a 
Significéent disadvantace. EFilimanation of this disadvantage 


N) 
Nt 





ie 


De 
that tke 


Wee-Or i 


< 


elemert vy 
Sequentiel 
weceecr is 


user verif 


It 
Teans cf a 
ClLecee Het 
Vaiue Entr 
JS) teh) eeatere 

e2 
wie. = iru 


variarle f 
an arkitra 
rumbker of 
Screen is 
thar 


liwe.teaticn 


=x; 


Cy. 


+) 


~~ mem 


*eractive 


progsen 
réequirec¢ input matrix or 


Sus CWLIineS were added to the 


ese matrix input subreutines were written such 


Peeceas flrst tntormed which specific matrix oz 


required; then premreted for the individual matrix 
These 
ly entered from the terminal. 
Peed 5 


mecrerOotTeance COrrecticn 


alues. values are then individually end 


@ree thesmattiy “of 
the terminal screen fer 


mt 2s returned =o 


if necessary. 
teractive sequential 
CC 


mane anc element 


data €ntry was pregrarr<eda by 


<wo-Gimensicnal TOCS Wath a <Grlanel promet 


eX 


Pesie2 on Dracr to the elénen. 
READ 


= 


Data element input vie a free-format 


Tem statenert. 
ce the matrix data entry sequence is ccmplets, 

ai igl 
with 


the maxinun 


$= returned to the terminal ecreer 


= we trix 
for ne SOW Ldem err lca TiCn.« 
myecaca fLicle@wadth ct 12 


fieteo x elerents avaitablie 


CrMat HWSecL Ease 


characters, 


On an 8Q column terminal 


Six. ERO Geue- Homme ce a ecolunimwGd. mension 22S Less 


this restriction presented no programming fcrmar 


— 
— @ 


tad 


NN) 





Basi. C2. Ras Cee aisLewe. ch CaS Wee see 1S a 
mcteas Le “OUid expect to sez it, providing the advazzags 
mom] 1n LOW ard cclumn identification. 

Nemes, | CPT SYSKy Subset ge iaatAl NL Sten eemas  ) ohas 


yvariatle-icrma 


ct 
Lf a 
rr 
=) 
ts 
ct 
“ 
(hy 
p62 
6a 
‘ny 
ty 
o) 
iD 
a 
t-} 
ae} 
(LD 
Lr @) 
ie 


Seweccarily «e2rninated wii Ck eoas; 


presuming tkat us 
é 


select alternat 
eee eG eat ezactive Input 


ies Wie = Lon s- 

Ae ohare tc reiieve the us2r of the burden of 
repeated system entry in th? interactive noda, several addi- 
tional gregram fiow ccntrol flags were adaed, allowing the 
GeeeereG: Saverng the entire oOriginai system of matrices for 
subseguéent calculaticns. Separate opticns for saving ¢ach 
SySt¢m@ matrix are attcmatically crfered at the end of ¢€ach 


beci@ran run. 

Tkese matrix saving options provide a further acven- 
Sagemro the USez in that the matrices are redisrlayed icr 
Sememeecaticn Dricr t¢ calculaticn execution. Users may then 
change incividual matrix elemerts, relieved of the burden of 


Mees yotem Te-entry. 


5- User-Defined Input Filss 





Perereughn thie kasic ockEjective of this werk is ‘to 
rrovide the user with a totally interactive method cf data 


input, séverai disacvantages tc the methcd of individual 





a * ~ - 3 
eee che Tent input are éprarent--input cf very tAarce 
_ —_- «= = = ° ~ - ad os ~*~ -_ * « mea re 
ieee ts aw cld¢y and time-consuning; inpuc cf systems o2 
-— - ~ et a ~~ ~ = $ 
feeeweceo WhCEe ECLSMSltS remain urcnanged from run <c xcun ie 


iieecmdcme tC FrOVide€ an increase 


aq measure cr user 
Pexabaiaty in data input, subroutine SETUP was modified to 
Peete werOvISiOrs f¢r matrix data input from @ date file cn 
@iecmcesa*S Gisk. Pieee 222s Shee Nac haces Wet Fae G cena 
ome a | May 2BOw bes iAnput [rem the user's disk. gal, sole pe 
PeocreauenCGuf2cetion is required of ths user as follews: 


u 
@., FRICMS Filedaf commands must be modified cr 
Bac oo Stew est 1c 


eacrk data sst. 


Ge the pane anda lWecaver. ce 


rk. The REAL Format statement (or statements) mtst 
Peechanged te rcretlece ehe pzcpst data fdoriac 


Greener User’ Ss Ineut file. 


Memcer tc en&> Tsasure ci flexibidztcy, che documenta- 
Efee Wheban Ssubrcutine SETUP was expanded to include several 
Specific examples of internal matrix data generaticn. The 


tnreée systenr matrices [Fj, [Gi, and [GAMMA] may generated 
either within user-sritten tweedimensionel DO-locps or by 
direct assignment statements. These methods may be prefer- 
able for the input of very large matrices with few nen-zero 
elemerts. 

Bh specific exemple cf internal program generaticn of 
SiemeGterte equation [H}] matrix is aincluded in sukreuvine 
BEALE. This matrix input methed may be preferable for the 
entry cf a large output equation matrix with very few nen- 
zerc elerents. 

Cnceé these acdificaticns nave been made? +t0 subreu- 
See mehttr Or Subrottine READE (as desized), the pregram 

nn 


snhoulé te re-compilec and then run in the usuai na 


Le) 
cn 





Mie icaeive PELOGLam trempt is provided at the beginnizg of 
Peewee ctfering the user che option o£ speciiying che 


sa 
Sipe. “CSCW cc GL tee =f eee. Sees iS 
BemereeO che {| Hj matrix (Or ¢cther Ffequired nous mateices) 
from s¢parate data files. USEIS With cudimentary prcaraa- 
Beg Skills may noweredify sukroutine READH (ce cme cf the 
Set eeecoecitie REAL subroutines) in the maimrer previcusly 
@eseraiced fcr Subrcutine SETUP Or subroutine Rea eH . 
Detailed examples of the nature and extent cf thase mecdifi- 


Seeens Rayebs fcund in Appendix A. 


7. Yate Extry Ccrrecticn 





imeem eerroOrt sto pretect users from er=>ors in data 
ap Ue, em G=rOr Correction sequence was incorperated into 
€ach wétrix input subroutine. 


Cnce the entire matrix cr vector is disp 


() 


a 
termiral screen the user is prompted with the q 
Vem eich tc change thre value cf any matrix élsn 
eereerue Ne*." TS the user types "No", progzta 
contirues. . 

iimeae User types “Yes", he is then proitfcted with 
three additional statements specifying the row number of the 
e€lemert tc ke changec, tne cclumn number of the element to 


MeO thet Mac tix 


Bowecenecd, “and the value to be inserted i 
element. After the corrected vaiue is entered, the new 
matrix values are returned Ee. the screen [o> 
re-verification. 

This correction sequénce continues indefinitely 


until tke user signifies that no additional changes are 
recessary. Program €xecution then proceeds nermally. 


ZG 





SY 2 = aS ~—- * - “~ ww —_— ie oe a cee _ et ~ 

nan ¥ te eeecclve CUOLDUAe Er ProgiaNds Secs. eae de G 

a a eee ee — we > Ee ~~ & —_ “~~ * em ow q ss om = = = os, - -~ = » eS —_ 

ae le ce ne ie ee ww tee Shae 2 re oe Gi eee ie aera nanatis 4a { ¢ 22 ww wae 

~ ~ » ~— ; _ —- a ~ ~ _-_ —_ bred 

Becouve rt Vv!) ShOULIag the user Inadvartenctliy make @n errenscus 

~ ~~ = = ae Tai ~ at = oy a ~ —_— a = —-_ ee om = 

kKeycoérc entry. PAcupee = er = Ness Wom eommen. oot = Crs 
iscclucde--Eentr ie k TRAa +A 3-5 ~+e an fom) Spe kara a a = _—< 

pe Se SG ww 36 cy e) a Bee en Che 2 Cc SIS Nad om Crewe G - & _— BOS Oe 

o eS Fe ~~ ww = on = —_ 

woaen tt! numerical vaius; EMenv. c 2 


tection features were 


Three sukroutines--KDREAL, RDINT, RDCHAR--were edced 


Pom@memoroX Jag esder to ensure that the froser input cete 
Measceeric Erevided tc the frogtan. SuUbEOucT oer URE sl were 
MeemmetG@icuwec:Y PCliny <€ real number or zere integer inrutnsay 


ke erccurtsred; subrcutine &ELINT is called at any feint a 
non-zerce integer input is required; subroutine RUECHAR is 
Supeegmet any POIlNt a@ Character string ("Yes" or "No") input 
Peomeeguared. 

Kithin each cf these subroutines a null string entry 
Protcecticn loop was inccrpcrated (allowing cne reccvery); 
Peometing the user fcr the correct data type input, and 
returninc an error message in tne event an incorrect dat 
fwocmas erecuntered. 

Within svbrottine HDINT, improper data “tyre entry 
N@emptieccke: OC rectec by the addition of a three-way IF 
ccmpariscn ci entry integer magnitude. Mavs MOditiceclon 
precluces illegal (bet automatic, with some ccempilerse) data 


Eve] GEEversSion ¢rrors. 


N) 
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iicseevecgrat design Ieatures boast the eiasmicre- 
Femeemtace Cr alilowirqg normai procram <csérminaticn 2: any 
Petpet cme Gets inpucw erase Ey merely Sreesing £26 “Er7erc" 


ma Se Ovace uUssem SErcr messaces in the event cf ineccnsistent 
BeccpegeelcGweCCNtLC, €lags, but te@Meratted the fErecran. 
Mees features was unrcesirable from the standpoint ci smocth 
Peeerce-lve Program Creraticr, and was impreved in OFTSYSX. 
Sieseucine CEFCK was mcdified <tc include EETURN 


Petementc any time inconsistent program flow contrel fi 
are ercctintered. The User iS rotizisd ci the type ci er 

encountered; that run térmiration has occurred; and yprespte 
Pegatdinc his desire to return to the beginning of the 


rrogrem cr terminate ¢executicn ccmpletely. 


Ze 





IV. ERCGRAM USE AND EXAMPLES 


[> — —_— > DP <a ap ee ee 


reccrc¢ed termina 
eGceenezel users should exasins carefully the éexa 
rrogram failure fo NM SeGelOn Ds Hews  esxea nol i 
Gemonstratées that Smenecdecs © (sneceseet Cre. ce ‘a 
certain design parameters may cause program fail 
tmeenrce: CULpUt!), ¢éveEen with the hignily stable "Qa" algo- 
Tithe. It aiso indicates cne possible method of correcting 
this tyre ci failure Fy merely making a very small change to 


cne cf tke design parameters. 


Ae CEFEN-LCCP ELGENSYSTES ANALYSIS 


The ‘fclliowing cren-lcop e¢igenvalue example was taxen 
boom [ Ret. S$, p.-669 ]. 

Exeeeca.i0Gn Of tre follcwirg frogram cutrut shcws cpen- 
ieee 8c “myedUCS at -t,0.-2, and -3. Wote that the sigenvec- 
Semcmor the Lert ang right €igenvector matrices (fe. 23) 
Gemccerord in cclumn fashicn te the open-lcop eigenvaltes 
Saeurated Mhhediately above them (fa. 32). 

The full terminal sessicn is recorded below, with user 


imeem, “cer Case jetters folicwing sach 2?" 


recore cr 

BEGIN RECCKCING CF TEFMINAL SESSION 

Pemetaeted 170.02 21:49: 

feeder Co term (recfaefa biksize 133 

glokal +xtlib fortmoéz mcd2¢eeh imsldp neninsl 


Wead ceEtsysu (start 
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Boece IiCk SEGINS... 


SenSyex 2S A COMFLETELY INTERACTIVE OPTIMAL SYSTEMS CONTRLL 
FRCGSAN. IT WILL SCLYE XUMZRCUS CONTROL PBOBLZINS OU TEE 
FCLLCSING TYPES CE SYSTEMS CCNTROL FQUATIONS: 

XDOT = (F)#X + (G)¥*U + (GAM) * (+496) 


MEASUREMENT ECUAIION-- 
Z = (B)*X + (D) 44 
REGULATOE PERFORMANCE INDEX-- 
(Yo#(A)¥*Y + UL¥* (8) ¥0) DT 
STATE FEEDBACK GAIN DEFINITION-- 
C = ~(C)*Xx 
[CC YCU WISH TC CONTINUE? TYPE "YES" OR "NO", 


Cs 
i 
5 
“\ 
Ae) 
+ 
bt 
A 
2 
tri 
Q2 
ri) 
em 
tr 


yes 
eee Ae Ne fo 
MimpeeeloueOrISYSX [5S SPECIFICALLY DESIGNEL TO REAL 
weitere x LATA INTERACTIVELY, SEVERAL ALTERNATE 
repo Ds ARE AVAILABLE TC USERS: 
Pees tat oe MEN YG", 6 AND "GASMA" MATRICES 
ewe SE SPE AD EF ROY SEPARATE DATA FILES. 
aah econ e FMW E". AND "“"GAHUMA" MATRICES MAY Bz 
Peeper CaityY DEFINED WHETHIN SUBROUTINE "SETUP". 
Piette ELTHES CASE, T8E USER SHOULD CBTAIN A CCPY 
OF THE EFCGRAM LISTING AND fcXAMINE 
THE EXAMELES CONTAINED IN S/he ore Ee.) 
poem cum L Sn ©£OC CENTINUE? TYPE "YES" OR "NO", 


yes 
Pemye@urmron LTO JNEUT THE "YF", “Gt, AND "GAMMA" 
Pernecrs FERCN SUBROUTINE “SETUP" IAW THE 
METHOD DESCKIBED CN THE PREVIOUS SCREEN? 
Pepe Yes sO R TeNOr, 
rc 


GENERAL OFTSYSX OFTIONS: 
GPTr@ie 1 —-—- SYSTEM ANALYSIS WiTHOUT 

CFEN-LOCP EIGENSYSTEM CALCULATIONS. 
@eemene2 =—~ SYSTEM ANALYSIS WITH OFEN-LOOP 
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ss) 


(wv 


PUGswoc STEM CALCULATIONS: 
CEP iOMP Ss =-— CEFR-LOCE ELGENSYSTEM FOUND 
END ee Case a eRe TNA TES 
Geo een tree Phi eY FOLEOAS TRME rear 2) 
Sew tea = VCLAL DISTRIEUTION MATRICES COZZEUTEL 
ROUT heer Oh REGULATORS S Pht EESes 
Creoboa iat ets ANALYSIS. 
SeovGl AN Of Pre. t).2,5, Cx 4. 


Wives + OF STATES {NS) OF THe SYSTEM MATRIX 
wo atrarR TX . 


Penayzemearol ER SETTINGS FOR TdIS RUN ARE AS FCLICWS: 


1G) Gay 1G Rec oerr OLTFA, ITF2 IPRS fEDEN LE ibEEUG 
2 0 0) 0 C Q 0 0 0 0 0) 
Peer eiae GPSD iL¥U iNCRM IREG NS NC NOB NE 

9) Q C 0 0 3} 0 Q 0 
Oaweeecr SYSTEM = 32 
NUMEER CF CONTROLS = Q 
POM neGr SCBSERYATICNS = O 
MeMeeen CF EROCESS NCISE SOURCES = 9Q 


6a) 


ENIER THE SYSTEM MATRIX "F"-=MATRIX 
DMS TON = # STATES (NS) X # STATES (NS) 
TEE ELEMENT F( 1, 1)= 


Poe eELEMENT F( 1, 2)= 


Meee clenoNT F( 1, 3)= 
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(HE ELEMENT F( 2, 1)= 


C 
TEE ELEMENT F( Zz, 2) 
0 
THE ELEMENT F( 2, 3)= 
TEE ELEMENT F( 2, 1) 
-6 
TPC rEMeNT F( 2, 2)= 
-11 
TEE ELEMENT F( 2, 3)= 
-6 
THE SYSTEM MATRIX “F"-MATRIX ... 
Cac '  4.0CC00 0.0 
Gao 0.0 1.00000 
—emecgod  =11.9CC00 -6.C0000 
DO YCU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO". 
Nn 


Cine LOCE DYWANECS MATRIX. ns wew cece ncscce Fae 
0.0 Os 100CE+01 0.0 
0.C 0.0 0. 1000D+01 
puecoueer+O1 -O.710CE+02 -0.6000D+01 


CEP tOCE EIGENVALUES. cee eee -DET(SI-F).. 
mieouugurtoCs—2. 000 00L+00:;-3.00000D+00: 





SOrew oO AUER Tg gtielve do Die yee 
c 


pee 64 — 2. fe2 17SD-07 1.0482¢50-01 
Se) FSSCs0-Ci 8.364 3258D-01 -3.74438550-01 
eee Stl -Ut —Sa7256716D-01 $.4345640-61 


Mee mecr Left EIGENVECTCR MATHIX...-..T-INY.. 

poem + CO - 423301270400 —8.660254D~01 
T2374773D4+C1 1.833030D401 4.582576D+00 
9.535$3925+00 %.430S5C9D+01 4.769696D+900 


OMPLEITE. DO YOU wANT ANOTHER RUN? 
Tere a Or UNO 
LO 

meeererer bots A 15 NCW TERMINATED. so. c0 0 


remeene 7 670.89 2171:52:C8 
mecorda cif 
END ERECCFDIANG OF TERYINAL SESSICN 


Po SO EGUEATCR SYNTHESIS 


Tre fcllowing regrlator synthesis exanple was taken fren 
OWeeeuLe Ncx~es cn Acvanced Centrol Systems", by Prefessor 
DemiameeCGli2ns Of the Naval Festgraduate School, Monterey,Ca. 
This example inveivec determination of the optimal requlator 
Gains kased on an arkitrarily chosen quadratic index; with 
the varicus system arc cost matrices described belcw. 

Exanination ci the extensive program output indicates 
Pieemenerecttanal Cregquiator gains ares ~=-5.0 and -SQRT(1¢.C). 
imemebectra.:c Sign cf the gains is consistent with the defi- 
rizticr digsrlayed on the first screen of the program (fe. 
34). 

Mr= tuli te 


Giinal sessicn is recorded below, with user 
input in leser case i¢etisrs 


POMEOwL EG seem 


Meccrigq Cr 
EEGIK EECORLING CF TEFAINAL SESSION 


Se 





ie 
é 
4 


Merece= £6 =erm (recftr fa bikedze 133 
el 


Peeeeeiwen stl 1b Ecrtmocz mod2¢ech imsic 
e 


C202) 13; 55926 


c 
CII 
¥SxX IS A COMFLETELY INTERACTIVE OPTIMAL SYSTEMS CONTRCL 
FFCGERM. IT WILL SOLVE NUMEROUS CONTROL FROBSLEMS CN TEE 
FCLLCKING TYEES CF SYSTEMS CCNTROL FCQUATIONS: 
XDOT = (F) 4X + (G)*U + (GAM) * (ii+%0) 
MEASUFEMENT ECUATION-- 
Z = (H)*X + (D)¥*e + V 
REGULATCR PERFCRMANCE INDEX-- 
J = 1/2 * INTEGRAL (Yt*(A)¥*Y + Ut*(B) ¥U)DT 
STATE FEEDBACK GAIN DEFINITION-- 
U = -(C) *X 
DC YCU WISH TC CCNTINUE? TYPE "YES" OR "NOM, 


--TATAa ENTRY-- 
BRITECUGH OPISYSX IS SPECIFICALLY DESIGNED TO REAL 
ALL MATRIX CATA INTERACTIVELY, SEVERAL ALTERNATE 
“METECDS ARE AVAILABLE TO USERS: 
METHCD 1--IHE "E","G", AND "GAMMA" MATRICES 
MAY EE READ FRCM SEPARATE DATA FILES. 
METHCD 2--THE “F","G",) AND "GAMMA" MATRICES MAY BE 
EXPLICIILY DEFINED WITHIN SUBROUTINE "SETUE", 
(NCIE: IN EITHEF CASE, TEE USER SHOULD CBTAIN &A CCPY 
OF THE EROGRAM LISTING AND EXAMINE 
THE EXAMELES CCNTAINED IN S/R “SETUP".) 
[O YOU WISH IC CONTINUE? TYPE "YES" OR "NCH, 


EG@eeroUp GSH TOeINFUIT TBE "FN", “G", AND “"GASMA"™ 
Miner CES FECM SUBROUTINE “SETUP™ £AW THE 
MEEHOD DESGrIBED CN TH PREVIOUS SCREEN? 

ieeD 2 Mey SU Oe ENG < 
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CELLO N 


("PU-MATEIX ENTRY 
CETION 4 -- MODAL DISTR 
®ITHOUT FILT2E 
SELECI AN OPTION: 1,2, 
7 
Z 
CC YCU DESIRE RYE YALUES CF 
TYPE "YES" CR "NO", 
EC 
CETION 1 -- CETIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 
CETION 2 -- OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 
CETION 3 -- CPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 
CETION 4 -- OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 
Sameer AN CETION: 1, 2, 3, 
? 


GENERAL CFTS 


bo 
' 
j 


OR RE CURA ZCR Ss 


OPT Le S.< 


4IfeCcuT 


5 
BIGENS LS 2 Eee Oss . 


pha ioe are nies re Nee efor wae 
Wi ae Gr got EA kee 


CALCULATIONS. 
PEGS SS fa 
Vote iN aS eS. 

PODLONS Aceon ae 
LBUL LION 


S CUD 


i 
He Cen eet 2.8 
4 


On SL tap aol aaeANe LIS LS. 


Si {Oey a 


oTaTE AND CONTROL? 


Cris towmbo hy Cia ssltCAL GFT LONS 


ANCT/OR REGULATOR 


NORECLe a NAL CU LOR th" 
AND/OR REGULATOR 
PAN Ai Ct 


AND/OR REGULATOR 
Paco k al, K's 


AND/OR REGULATOR 


PPE ate E CANE SK" 


OR 4. 


fmemeueteewiSH +O LETERMINE THRE STEADY-STATE RESFONSE 
PeneA C@NStAaAhl DISTUREANCE? 


Shs 





Peek wits CREO NS 


rs 
© 


Te a te ee OAT OSS Ta Le ee 
aie “Ys MA dd oe 8 dee » aad wa Zs AL Des ie TEU ieee 


4 
om) 
ms 
©) 
cS 
+4 
4 | 
Z ats 
| 
7 ee) 
t-3 


ieee Ss Cnesot. 
LC 


Cra NeLOOL TRANS FeR FUNCTION OPTLONS: 


eee te => NC OFPEN~LCOCE DRANSPAR PUNCTIONS COMPGT ED. 
Wempen 2 =- PCLES, RESIDUES, AND ZEROS COHPUTED. 
eres =- OLY FrOLES AND ZEROS COMPUTED. 

Stuns == OLY POLES AND RESIDUES COMPUTED. 


See GTeaN CETECN: 1, 2, 3, ORs 


é4) 


Nom ot TRANSEER FUNCTION OPTIONS: 


Sweet =>—- KC NOTSE TEANSS ak FUNCTIONS COMNPUZEL. 
Seton 2 —-— FLEES, RESIDUES, AND ZEROS COMFEUTED. 
Greeeouess——- ChIY POLES AND ZekOS COMPUTED. 

CeT ION => Cui ee Ores eat esc ot DUES .CONRUTED-. 


4 
SeeeuCrVaNeCETIONS 1, 2, 3, Or 4. 


G@eePeNSATOR TRANSFER FUNCTICN OPTIONS: 


SereeeN 1 ==" hE COMP. TRANSFER FUNCTIONS COMPUTED. 
Coimmeime 2 == ECLES, KESITDULS, AND ZEROS COMEUIED. 
Saerowes ——= ChEY “POLES AND ZeRCS. COMPUTED. 

Sev TOveea—-- GLY POLES AND RgeSEDUES COMPUTED. 


(NOTE: A CCHUPENSATCER TRANSFER FUNCTION CAN EE 
COMEUTED CNLY IF BOTH A REGULATORAND 
Clb nan tro ent Ae SEZ EL Sai, Che INPUT.) 
SeeEeL AN CETION: I, 2, 3, O02 94. 


ea) 


MILL A FEED-FORK@RD DISTRIBUTION MATRIX 
web aM ATRIA) BE LNEURS? 
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Peep oe het yO. 


THE CRDER CF THE NUMERATOR ECLYNCHIAL OF EACH TRANSFER 
FUNCTICY. 
ALL "Ni ZEROS OF THIS PCLYNOMIAL ARE 
Memes PEoT TEILS EC MANY EXTRA RGOTS 27S 
HAN 10.0¥** (-I5) IS CCNSIDERED ZERO. 
EE LEFAULT VALUE CF THIS FARAM=TER (IE) IS 6. 
IN CTHER WORDS, IE = 1.CE-6. 
IF YCU DESIRE A CIFFERENT MARKOV CRITERIA, TYPE THE 
INTEGER VALUE. 
IF YCU CESIRE THE DEFAULT YALUE, TYPE "O" (ZERC} 


DO YCU LESIRE TC SYNTHESIZE A STABLE FILTER (OR REGULATOB) 
BY DESTAEIIIZING THE CRIGINAL SYSTEH? 
(NOTESWORKS FCR FILTER OR REGULATOR BUT NOT FCR EOTH 
IN THE SaME RUN.) 
TYPE "YES" CR "NO", 


Lo 
Pemreumwbhothes LC PRINT THE EULER-~LAGRANGE EIGENSYSTIER 
ERIeCR £0 DECCEFCSITICN (FOR CHECKING THE PROGRAM) ? 
Mae Ao! OR NO. 
yes 


Merecers-fCIRAL LENSITY (PSD) OPTION 7 ; 
mere tb ——- COMEUTE THE FSD OF THE OUTPUTS AND/CR THE 
CONT nor osOret as —CONTRODEEE SYSTEM AEEN 
POMmGED EY FROGCESS AND MEASUREMENT NCISE. 
Miome > #5Glad A REGULATOR ANE A FELT ER MUSTEs 
Roe Drh eeN THESEROGRAM  fO USE THIS CPIICh. { 
Grime n 2 =—- SAGE AS GCETICN 1 ABOVE BUT CNLY PRINT TEE 
wc ivUts «Gre rnecue rn ANor ER FUNCT ION 
Sten Eee oD COMPU TION. 


=) 





merch Ss == NOT CLESIRED. 
= Ve 


Wroe. AN CeETION: ee ron 


CC YCU DESTRe REGULATOR SYNTHESIS ONLY? 
ieee Yoo" CR UNO, 
moe rk Ge ae OF STATES (NS) OF THE SYSTEM MATRIX 


Gee PATRI). 


Meester ¢ OPNCONTEOES (NC) OF THE CONTROL SYSTEM #CIEL 
Goc =EATHIX). 


@a) 


ENTER THE # OF MEASUREMENTS CR OBSERVATIONS (NO) OF THE 
("BEN ATRIX) « 


ea) 


W) 


ENTER Itz #4 OF FROCESS NOISE SOURCES (NG) OF THE 
('GAMMA"-MATRIX). 


e,) 


PBC/EAPFANETER SETTINGS fOR THIS RUN ARE AS FOLLCAS: 
immeceman C55 LM IfFl ITF2 itF3 IFDFW IE ILDEBUG 


1 0 C Q C 0 Q 0 0) Q 1 
Meer DStAr L£PSD IY0 INCHY IREG NS NC NOB NEC 
Q C Q C 0 1 Z 1 Zz 0 
Cees cr SYSTEM = 2 
NUMEER CF CONTRCLS = 1 
 NUMEER CF CBSERVATICNS = 2 
Miers GEeEROCESS NCISE SOURCES = 0 


Meelea Pee SYSTEL MATRIX “F'~-MAT AIX 
DIMENSION = # STATES (NS) X # STATES (NS) 
MereerGrMeNT F( 1, 1)= 
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Te SlewontT F{ 1, 2) 


TEE ELEMENT F( 2, 1)= 


0 
eemeELEMENT F( Z2, 2)= 
: 
ice ome ee Ana ke Pm MART ae 
50 1.0CC00 
C.0 = 0 
DO YCUC WISH TC CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
eee YESH CR "NOY. 
LO 
Pearce OP YNAMLCS MATRIX. scccussccenesscclac 
o=¢ 0.10 00L+01 
0.0 0.0 
Ger LTCCP EIGENVALUES... ....-+«-DET(SI-F).. 
0.0 Barre V : 


Creel GCCP RIGHT EIGENVECTCE MATRIX. ...22501e05 0 
imu cuu@cl+O0 =-1.000CCCD+00 
0.0 2-«220446D-16 


Ghee LOCe LEFT EJGCENVECTIOR MATRIX...-...T-INV.. 


TOCCCCECD+CO 4.503 6C0Dt15 
0.0 4.503€COD+15 


Diener EASUFENENT SCALING MATRIX "“H"-MATRIZ . 
DIMENSION = # CBSERVATIONS (NO) X # STATES (NS) 
Beer BMENT H( 1, 1)= 


#0) 


a9 





THE ELEMENT H( 1, 2)= 


eect h iment Hl z, We 


6 
THE FLEMENT H( Z, 2)= 
: 
TEE MEASUREMENT SCALING MATRIX “H"-MATRIX ... 
Pecos C 0.0 
0.0 1,0C C00 
CO YCU WISH TO CHANGE THE VALUE OF ANY HATRIX ELEMENT? 
TYPE "YES" CR NO", 
ro 


eee Pee SCALING MATRIX 2. ccs sec ce eee ole w 
eee COLr+ 0 1 0.0 
0.0 O2100CI+01 


ENTER THE OUTPUT MEASUREMENT COST MATRIX “A'“-MATRIX . 
CIMENSION = # CESERVATICNS (NO) X # GESERVATICNS (KXO} 
TEE ELEMENT A( 1, 1)= 


z5 

THEE ELEMENT A( 1, 2)= 
C 

TEE ELEMENT A( Zz, 1)= 
0 

Teo eLAee NT At Z, 2)= 
C 


GoaeeeumecGl MEASUREMENT COST MATRIX "“"A™N-MATRIX ..«- 
Zee OU GC Oo & 


4G 





Be YCG WISH TO CEANGH THE VALUE OF ANY MATRIX ELISE: 


ooo ars CR NO, 


ee 
COUTEUT cost DR ete re ta er ce nc ko oe See hens 
C.Z5(0D+C2 0.0 
0.0 0.0 
Mien TH® CCNTRCL DISTEIBOTICN MATRIX “S"-4ATETX . 
DIMENSION = # STATES (NS) X # CONTRCLS (NC) 
TEE ELEMENT G( 1, 1)= 
0 
TRE ELEMENT G( Z, 1)= 
1 
TEE CONTROL DISTRIBUTION MATRIX "G"=MATRIX ... 
6 
sPROECLOLOR® 


Pemereumearch LTC CHANGE THE VALUs OF ANY MATRIX ELEMENT? 
ieee Yoo Chat NO™, 


LO 
EXTER IHE CCNTRCL COST WEIGHTING MATRIX “B"-MATRIYX 
DIMENSION = # CONTECIS (NC) X # CCNTROLS (NC) 
TEZ ELEMENT B( 1, 1)= 
MEE MGONTEC 1 COST MANET X sans --cs.s.000 Bo: e's 
1.00000 
DO YCU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 
ro 


TEE CCNTROL LISTHRIBUTION WNATRIX.. 2222.2 oe -Gee 
0.6 


O. 10CO0L+0 1 
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fereerwimol COSf MATRIX .c<ccecceeescccecneeBe s 
Se GC COLtS 3% 


Bieta aANGh SYSTEM MATRIX... 


0.0 WeOO0CCCODECO O20 0.0 

Ore 0.0 0.0 —WacwouOGCr +00 
eaeoueuccD+ct 0.0 ae ee, 

0.0 0.¢ -1.C00C00D+09 0.0 


Gero Of Ahe Ziok. 24 


PrGenvelLURS 24ND EIGENVEC' 
pe as HO2 2.028 s¢eu@eee eee 


PULERroRAGKRANGS SYSTEM 


pees sD+CO 1.587139D+00 
mere sch +GO= 1.59112S5D+00 
ees che CO 1.981139D+00 
oro |ocpeeO- 1.5681125SD+00 


ceeeeowtah—G2 7.168812D-02 -1.824925D-01 -1.503482D-01 
Preeti eo 2e508245D-01 -5.082459D-C2 =-5.262675D-01 
mem eee —-2.0683770-02 2.631337Dt00 -2.541225D-01 
eereoreeecr =O! —~5./15222D-01 -7.517412D-01 9. 124627D-01 


PiiGeonetctaad CF CETIMNAL REGULATOR. 2.2. cease 


POPC PY ECTORS FRCY RGAIN FRIOR TO CNORM 
oememes c22)-Ci —2 .6 16312D-03 
eee izb~O! -2.26657./D-01 


GeveGrecrlti Mal REG. Em~VALUES. .. DET(SI-FtG*C).. 


C-LCCE BIGHT EIGENVECTOR MATRIXec eee ee eee Meee 
womete ze cp! —3.162278D-01 
moCGCOGeb!eCO 0.0 


MGMrrcr EI GENVECICR MATRIX. .cccses occ e e CK. 
—lIeccutsertO0  1.2€1139D+00 





Meeeeer cel. GEG. LEFT E-VECTOR MATRIX..sM-INV.. 
Or eo Ue GCC OFC 
1 


eee eet —-1.000CCCD+00 


qs NEU Ale ee PS or = Oe tee (ee se re Win) Oo) Ge ee Ne AS nee 
Ri ¢cedirsatK Gain Mees, itwerctc Ae a 6k. =e es. © Se 8 a 


Boece lOsey LCOP DYNAMICS MATRIX ...-0.F°G¥C.. 
0.0 T.000CCCD+990 
pea JCCUCCDECO —-3.162278D+00 


Pobsti cis COMPLETE. ©O YOU WANT ANOTHER RUN? 
ero Seuss Oh “Net, 
no 
Sitter ooo kok TS NOW TERS LNATED. «00 ces 


peer =ol4271.85 14°03:C3 
MeeGro off 
END BECCECING OF TERMINAL SESSIO} 


ee fer eR SYNTHESIS 


Tke following Kalman filter synth2sis example was taken 
fren NLecttrs Notes on Advanced Control Systems", bY 
PeOrcccor [.J. Cellins of the Naval Postgraduate Schocl, 
Monterey,Ca. 

This example invelved determination of the optigzal 
filter gains of an arkitréry system; modeled nearly identi- 
Swerve uGlethe previous regulatcr problen. 


In its present ccnfiguraticn, OPTSYSX program sequencing 


reguires the design of an optimal reguiator, shgmee ss eS) 
perfcrming ény ocptinzael estimatcr synthesis. Pie —Cnden. eo 
CGmply with built-ir program sequencing conventions, and 


Circumvent crogram difficulties which may not be specified 
in the particuiar system model, optimal filter synthesis may 


re acccrplished Ly ertéering the identity matrix [I] in these 


fox 
tar 





wer 


(y 
ra 
5 *o 


deed Ss =: @o 
ee a hed E 
ay 


EG y 
aects ali pe Cesrec> Por tie Syor- 
ext 





f 
+ 


mcdei. 
€nsive program output indicates 


C 
na 

Siac the Orctimal Cillter gains are: -5.0, and -SQRT(2.0). 
ul 





Tre full terminal sessicn is recorded below, with user 
Mmpuc an LoWer case Jetzers fclicwing zack "2?" . 
mM=CCre CF. 


c ; 
BocitMPbECOALING CF TEEFAINAL SESSION 
Pee i=C.0170.02 21:493: 

Beeecer C6 term (rectim fa biksize 133 

Global txtilackb fortmoac2 ncd2e€éh imsidp ncnimsl 
dead ¢cprtecysx (start 
BLSCCIECK BEGINS... 

Oem oro A VCOMOLETELY INTERACTIVE OPTIMAL SYSTEMS CCHKTECL 
Pacer. be WLiT SCLVE NUCEROUS CONTROL EROBLENS CN THE 
MeMmeewliicetiPsS CF SYSTEMS CCNTROL EQUATIONS: 

XDOT = (F)#X + (G)#U + (GAM) *(a4+¥0) 
Peacunc icon EC ORTLON=— 
Ze= {jek + SD} ew + 
REGULATOR FPERFCRMANCE INDEX-—- 
Ue=si7z * INTEGHAL. (Yt *(A) ey + US* (8B) *U) ET 
Srp ere Ey BACK GALNS DEF TINITLCN== 
Uie= = (C)ixk 
ROMY OU AISt: TO CONTINUE? TYPE "YES" OR "NC. 
yes 
= eA A, Se cee 
PenewGnmeOrisySx 15 SPECTFICALLY DESIGNED TO REAL 
Mere atark LATA INTERACTIVELY, SEVERAL ALTERNATE 
METHODS ARE AVAIIABLE TC USES&S: 
METHCD 1--THE "F","G", AND "GAMMA" MATRICES 
MAY BE BEAD FRCM SEFARATE DATA FILES. 
METECD 2--THE "F","G", AND “GAMMA" MATRICES MAY EE 
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LO 


ea) 


no 


[CO 


Peele Crile Pak INEE atlas SUSROUTIN EW YSETUE" . 
con Seto yee Stk SHOU Opr Adelina Cord 


CF THE PROGRAY LISTING AND EXAMINE 


YOU WISH TO INPUT THE "F", “GY, AND “SGANHA" 
MATRICES FFCM SUBROUTINE "SETUE" IAW THE 
METHOD DESCFIBED CN THE PREVIOUS SCREEN? 
TYPE "YES" OR "NO", 
GENERAL CPISYSX OPTIONS: 
CETION 1 -- SYSTEM ANALYSIS WITHOUT 
CEEN-LOCP FIGENSYSTEM CALCULATIONS. 
CFIION 2 -- SYSTEM ANALYSIS WITH OFEN-LOOP 
FIGENSY STEM CALCULATIONS. 
CETION 3 -- CEEN-LOCE FIGENSYSTEM FOUND 
AND PROGRAM TERMINATES. 
("FM-NYATEIX ENTEY FOLLOWS IMMEDIATELY.) 
CETION 4 -- MCLCAL DISTRIBUTION MATRICES COMEUTEDL 
WITHOUT FILTER OR REGULATOR SYNTHESIS 
CR STEADY-STATE ANALYSIS. 
SeUEG@T AN OPTICN: 1,2,3, OR 4. 


MeUmBrotkRE RMS VALUES OF STATE AND CONTROL? 
myer so" CeO". 


CEES iste Pen Chasstean. CeLTiICNns: 
Sereoy t == OFTEMAL FILTER AND7JOR REGULATOR 
SYNTHESIS WittieNOmekLERNAL "Ci! OR K's 
MATRIX INPUT. 
CETLON 2 -- CFTIMAL FILTER ANCT/OR REGULATOR 
SYNCH boissons EXDTERNAL “C" 
MATRIX INPUT. 


4s 





ed 


pepe 3 => CP@eUAL ELLIER AND/OR REGULATOR 
SMT A ESTS SNES R Se ERNAL A" 
WaTaIlX INPUT. 

Ceaieeen 4 == CETENAL FILLIES Ane /OR MB EGULATOS 
Soe ei SAL ENN ee SANT MK? 
Baelki xX INPUT. 


ool AmCeL tO cei 2° 3, OR 4. 


femowesaeod TO LETERSIT NE THE STEADY-STATE RESPONSE 
FCR A CONSTANT DISTUREANCE? 
ieee is OCR NOW. 


Memerouewe so TO CETERNE NE THE MODAL DISTRIBUTION 
mele GaN ei AtTricks? 


eee Saks Ck MNO, 


Wee LCOP PT PANSEER FUNCTICN CPTIONS: 


Greet == NO CPEN-LOOP TRANSFER FUNCTIONS COREUIEC. 
ese s=—- POLES, RESIEUES, AND ZENRGS COMPUTELC. 
Geiaeh 5 == ONLY POLES AND ZEROS COMPUTED. 

Galen + == ONIY PCLES ANU RESIDUES COMPUTED. 


SELEGl AN CETION;: 3, 2, 3, OR 4. 


Neer TRANSFER FENCTION CETICNS: 


GEmmenve l ~-sNO NOISE TRANSFER FUNCTIONS COMFUTETL. 
Ore ON 2 —~= POLES, KESIDUES, AND ZERCS CCNEUTED. 
ries ies =—- ONLY EOLERS AND ZEROS COUPUTED. 

PePrON ys == CNEY FCLES AND RESIDUES COUPUTED. 


Speer AN CETIONs 1, 2, 3, OR 4. 


Gur owenrOn LRANSFER FUNCTICN OPTIONS: 
CETION 1 -~- NC COMP. TRANSFER FUNCTIONS COMEFUTET. 
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GET GN 2 == CLES, EESIDUES, AND ZEROS COMEUTECL. 


eee 5 -—— CNECY POLES AND ZEROS COMPUTED. 
See = CHhLY POLES AND RESIDUES COMDUTSED. 
(ereee A CCOAPPAS ATOR TRBiSscek FUNCTION FS002 


Cerne ChMieuei: EOhwy A REG ULe To: 
Be cee ae Sy NT RST Zep 
Bo ACh ei ru.) 

Soeseur ace LLONS Ty 2, 3, OR. 4. 


Mitte Foe D-rOnklhD DISTRISUTION MATRIX 
ees HATRIX EE INPUT ? 


Pipe reo! OR SNC tt. 
ro 
Perce eereeeO mh LoORMEbES THE CRITERIA FOR CLECIDING WEEN 
A M@RACV FARAMETER I5 ZERC-TRE MARKOV PARAMETER 
iivtoiermees THE CabDEr OF THE NUMERATOR POLYNOMIAL OF EACH 
TRANSEER FUNCTION. 
Peewee aOs CF THIS FOLYNCHIAL ARE EFRINTED OUT AND 
iieomierot TELLS HOW MANY EXTRA ROOTS EXIST AT Z = Q. 
eran 10. O° (-TE) IS CCNSIDERED ZERO. 
ipnee EA ULT VALUE CF THIS PARAMETER (IE) IS 6. 
fee intceewonvS, LE = 1.0E-6. 
memiounuscrneE A DIFFERENT MABKCV CRITERIA, 
ieetere’ INTEGER VALUE. 
feu oectan THE LEFAULT VALUE, TYPe "0" (ZERO) 
g 
C 
Pomeummosirne LOC SYNTHESIZE S2 STABLE FILTER (OR REGULATION) 
Diet eo cweret iis NG THE CRIGITNAL SYSTEM? 
OROTESWCRKS FOR FILTER CR REGULATOR BUT NOT FOR ECTS 
IN THE SAME KUN.) 
eo reo" CR “NO. 


LO 
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he, 


ea) 


(a) 


no 


ey) 


N) 


ea) 


EC yc DESTISE 


Ene Cs Te 


ies YES" CR 


Coren 2 -— 


Crierem 3 —— 
SoLoer 


PeexyeU DESIAS 


Pere 


ENTER THE # 
erV-WrTRIX 


Pern IEE # 
eo am ATR EA 


PrrpR Ihe t 
WH =m AaR Tx 


Pein Lhe ¢ 


"YES" CR 


Be SERENE THE £ULen-CAGRANGE. 2 2G 1S 1 se 
Deret(eOstreo le (FP Ox CHECALUG THE 


HNO", 


L Ce OO ee 
SCOUWPULE Ui eesD Or THs OUTPUTS Bil /Ge 
Pic eCONTRORS TOF THr CONTROLLED SYSTEC 
WoeN PORCEL EY FROCESS AnD MEASUREMENT 
nO oul (Horas BU Lit oa AWD a 
Seneneerwo | Sereno LDENT DN PEE SPROGEAS 
POG fits GE LION .} 
Sereno Ie PECnat ABOVee bY ONLY Phen TEE 
Root eonoeOr EACH TRANSFER FURCTION 
Uso cst Neth eee Sl CONPuUraArTai ON. 
N@m ee LESTRSL. 
ANGEL LEN? 1302, 


hat 


Sie Tal ohn Si a Cys eae 
x Lu pi Luar 


Ong She 


RecummOUGR SYilanets ONE Ye 
nVoOM, 


Orci aves — (NS) Of THE SYS PEN MATRIX 


Orme bole (NG) ORs tae SYSLEN MODEL 


OF MEASUREMENTS OR OBSERVATIONS (NO) 


Gee seGeos NGlst SOURCES (NG) 


MGM Me —AATRIX . 
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Meoyemrarsce tk SETTINGS FCR THIS RUN ARE AS FOLLCEeS 
MP~mneeice ios di Jiri iff2 ifFs IfDOFW Is Licus 
Q ‘ 0 C 0 0 Q Q C g 
eee see 6 P SD hUTYU6CUINORHS)6UETREG 6 6UNS)6UNC PROB UONG 
C Q C 0 #) 2 1 1 Q 
O@eeer OF SYSTEM = «¢ 
fem CF CCRTacLs = 1 
niger CPeCRSERVWATICHS = 1 
Peteeenmereer ROCESS NCISE SOORC#S = 9 


Peet heeotOLEN MATRIX “FW=MATRIXZ 


Pema noPON = « STATES (NS) X * STRTES (NS) 
Meret bheientT F( i, 1) 


oa) 


0 
THE ELEMENT F( 1, 2)= 
1 
TEE ELEMENT F( 2, 1)= 
0 
Tec ePEeRE NT F( 2, 2)= 
C 
TRE SYSTEM MATRIX "F"-MATRIX ... 
G10 1.00000 
bes 0.0 


DO, YCUt Wel SH TO CHANGE THE VALUE OF ANY MATRIZ SLEMENT? 
ieee eS! CR MNO. 


Lo 


CEZN LCCP DYNAMICS MATRIXw wc cee wee ecccece oF. 
0.0 elo OCrEt0 4 


4s 





5 
TEE FLEMENT H( 1, 2)= 
C 
TEE MEASUREMENT SCALING MATAIX "H"-MATRIX ... 
aererc occ 00 
CC YCU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEKENI? 
TYPE "YES" CH "NO", 
no 


ieee meme, SCALTIG*®BAP EI X. cect ec ewes esaatas 
0. 10 COT +0 1 Q.0 


Peer THE OUTPUT HEASURSHENT COST MATRIX “AN-MATRIXZ 


DIMENSION = # CESERVATICNS (NO) X # OESERVATICNS (NO) 
TEE ELEMENT A( 1, 1)= 


1 
eeeeeiee “CRA SUREMENT CCST PMWALTRIX “ANU =UATRIA ... 
1.€C00C 
PGweecU $aSH TC CHANGE THE VALUE Or ANY MATRIX ELEMENT? 
TYPs @rss" Ca NO", 
Zo 
CCTECT COST FATH TX. cee we ccc cece ec cc ccc cc ccs Ane 
0. 10COL+01 
Mmeene He CONTRCL  DISTRIEOUTION MATRIX "“'GY-MNATRIX . 
Pine wolON = 4PSTATES (NS) KX # CONTRCLS ({NC) 
PeeeeecioNt Gi 1, 1)= 
0 


3 C 





ra 
Mil 
tri 
tr 
t4 
nw 
tz} 
mr 
+3 
GQ) 
ae 4 
NY 
- 
ae | 
eeu 
Nf 


cS 


SoeeeerrnO: DasteieuilOn sare ix “G'=$MATRIX ... 


aU 
0.0 
Demet diol TC CHANGE THE VALUE OF ANY MaTRIkX ELEMENT? 
Peete! Cr "NO. 
Yes 
eer THE RCWeNURZER OF TEE BLOWENT TO BE CHAUNGEC. 
Z 
Beenie CCLUMS NUWBER CF THE ELESNENT TO BE CHANGEL. 
E 
1 
Toe ELV EMENT G( <, i)= 
Ei 
1 
PimeewenOl DESTRDsuUGtCN MALRIX, “GM=MATRIX .... 
C.0 
1. CCO00 


Demeeetewiodt Te CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
iY PE owen. NOM. 


LO 
(ene the CCNTRCL COST WEIGHTING MATRIX "BYW-MNATRITX 
DevosiON = # CONITECIS NC X # CCNTROLS NC 
toe Peewee NT B( 1, t)= 
1 


oa) 


TEE COATRCI COST FATEIX...--.....58..-. 
VeeG00 0 
Demveumngon (TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
fpr Ch MNO". 
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Pee ming DISTRIBUTION MATRIixX.-cccceccscGes 
O.C 
Cee rouge rry | 


Meeeoecu rol CUS2 fATRL A scccc as « rer ee svete) Dies 
Oe 1eCGr+0 1 


Peeve ost CF CEIITMAL FEGULATOR. .202s0<0- 


eee ora aN AL REC. E-MALUGS. «. DET(SI=-F+G6*C) .. 
meemueetcrr=O1, /.071075-01: 


meer ewan BIGEAVECTOR TOT hlX 2.0 .s cee ccs Mes ee 
—mommnooen—O i —7.071CEED-01 
lbeeecccCenpt+CO 0.0 


Oeireecr EIGENVECTICR MATRIX. snes cvesscc cel *He « 
seme vceeD-01 7.C/71068D-01 


eee, Chl. REG. LEFT E-<VECTOR MATRI£Z..M-INV.. 
0.0 1.000CCOD+00 
-1.4142740+CO -1.000CCCD+00 


THE CETIMAL FEELBACK GAIN CCNTROL MATRIX...C=BINV*GI*S... 
~TeOUGGEftOO -1.474215+00 


reeleshD LeSr LCYNARITCS MATRIX ......F-G*C.. 
0.0 1.000CCCD+00 
“1.0CCCCCD+C0 -1.414214D+00 


ANALYSIS COMPLETE. TO YOU WANT ANOTHER RUN? 
TYPE "YES" OR "NC", 
no 
cecces eeOPTSYSX IS NOW TERMINATED......-. 


Pome 2.64 15:36:45 
mecerca cif 
END BECCELING OF TERMINAL SESSION 


un 
ND 





CT. EZAMELE OF PROGRAK FAILGUEE 


1. = ~—. 


iMeueras CcXamp~le, the choice of the quadratic in 
ree ACCOR SCrOmCting ErcGtaga instability; heoaas 2g EC 
e C J 


Poreetud PECGLaM Laliavze in subroutine dQrxZz. fre 
3) 


BeQernec> GCceinS ck -5.1, and -3.1 (pa. 6 





Meee (sliiaghz:’ icditicaticn of the ccst na ne 
Beevacuc value cf 4.0000 toa new value of 4.0001, the 
BaoGawce We= TUR a SEeccnd tite. Failure did net occur cn the 
seccnd run, anc the new calculations (pa. Soy indicarcs 
filter gains of -2.0, and @ "small" residue of 3.19D-14 
(essentially z¢erc). These are the correct values. 

This example Beams Beni one possible methed cf 
emma cwarg Gar tein ercgran failure modes, shcuid they ceccur 
@izince executicn. 

Tre full terminal sessicn is recorded below, with user 


S 
mppie dn lewer cese letters fcllcwing each "2?" 





Hewtemng che brcgram failure exampie, that perticn cf 
the repeated terminal output was delated up to the roint 


where prcarém executicn of the second run begins. 


BeCCrd cn 

BeGlnerPCGBLING CF TEFYINAL SESSION 

Pee 170.02 2134935 

filedef C6 term (recfm fa biksize 133 

gQlotal txtlib forztmoé2 med2ee¢h inmsldp neniasl 
Héad CEtsysz (start 


Peeeu re kh “EEGINS... 


Opens ATCOMPLEZTELY INTERACTIVE OPTIMAL SYSTEMS CCHKTECL 
Ppe@GraAt. LE wLLE SCLVE NUHERCUS CONTROL PROBLEMS ON TEE 
PeEeonwNG TYEES OF SYSTEMS CONTROL EQUATIONS: 


aS, 





XDOT = (F)¥*X # (G)*U + (GAM) *(W+40) 
MEASUREMENT ECUATION-- 
Z = (H)¥X + (D)eW + 9 
REGULATOR PERFCRAANCE INDEX-- 
1/2 * INTEGRAL (¥ * (A) *Y + UG *(By*#Up OT 
STATE FEEDBACK GAIN DEFINITION-- 
U = -(C) *X 


Teaco 2omecnc ths? TYPE YTS" CR NC, 


q 
il 


bud 
() 
wy) 


Sab a baa ONG 
Poeeeemcumort SY Sk [5S SPECIFICALLY DESIGNED TO SEAD 
ALL MATRIX DATA INTERACTIVELY, SEVERAL ALTERNATE 
Peeeee Ls BRE “AVAILABLE TO GUSERS: 
Deepa = -THE “PN wet, AND "GAMMA" MATRICES 
ma SE READ FROM SEPARATS DATA FILES. 
Paemene 2 = -The WEG AND “GAMMA MATRICES MAY EE 
Eror LiCLILY DEF Dit OeemETHIN SUBROUTINe "“SETUE™. 
Wem N DLIEGR CASE, THE USER SHOULD OBTAIN A COFY 
Ot THE EXGGRAN LISDING AND EXAMINE 
fee cA PL EDS CCNTATNED IN@ey rn “Ser U Pe) 
Rem Omen iwoH TGC "CONTINUES TYPE “YesIeiCR NGI. 


bemrveumemrcd TO TNPUT THE "PF", "G™, AND “GAMMA" 
Werner s FFCN SUB FOUTIN:e "“SsTUP" LAW THE 
PerdOD DESCRIBED CN THE PREVIOUS SCREEN? 
Pee IY eS SO ne NOM 
ho 
GONE Mole OrTS YS’ OPETONS: 


SEeerete id == SYSTEM ANALYSTS WITHOUT 
a CEBN-LOCE EIGENSYSTEM CALCULATIONS. 
Ghia 2 ——= SYSTEM ANALYSIS WITH OPEN-LOOP 
PIGENSYSlTEM VCACLCULATZONS- 
Sabwoueso —~- CFEN-LOCP EITGENSYSTEN FOUND 


AND FROGRAM TERMINATES. 
("F"=}MATRIX ENTRY FOLLOWS IMMEDIATELY.) 
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weer ONn 4 == CL 
PH OU eee Roos REGULATOR SYN tH=aS ES 
Como teNt -STALTE. ANALYS@S. 

2 


See fot AN Olt ON: 1,2,3, OR 4- 


Pere tevnoths xrt€S VALUES CF STATE AND CONTROL? 
meee cS CR UNO, 


yes 
CETSYSX LOR/CLASSICAL OPTIONS: 
CETION 1 -- CETIMAL FILTER AND/OR REGULATOR 
SYNTHESIS WITH NO EXTERNAL "C" OR KM 
MATRIX INEUT. 
CETION 2 -- OPTIMAL FILTER AND/OR REGULATOR 
SYNTHESIS WITH EXTERNAL "ct 
MATRIX INPUT. 
CETION 3 -- CETIMAL FILTER AND/OR REGULATOR 
SYNTHESIS WITH EXTERNAL "x" 
MATRIX INPUT. 
CETION 4 -- OPTIMAL FILTER AND/OR REGULATOR 
| SYNTHESIS WITH EXTERNAL "C" ANE Kn 
MATRIX INPUT. 
SELECT AN CETION: 1, 2, 3, OR 4. 
? 
1 
[O YOU WISH TO LETERMINE IHE STEADY-STATE RESEONSE 
FCR A CONSTANT DISTUREBNCE? 
TYPE "YES" OR "NOt, 
io 
[IC YCU WISH TO LETERMINE IHE MODAL DISTRIBUTICN 
AND GAIN MATEICES? 
TYPE "YES" CR "NOM. 
ZO : 


Ceri LCOr TRANSEER FUNCITICN OPTIONS; 
CETICN 1 ~-= NO CPEN-LCOE TRANSFER FUNCTIONS COMFUTEL. 


SiS 





ete SteePOLES, SRESIDGES, AND ZSROS COMPUTED. 
eee) 3 => ONLY FOLES AND ZEROS COMPUTED. 
Wepece 2 == ONLY POLES &A 


N) 


Greet TRANSEER FUNCTION OPTIONS: 

Seem -—- NC NODS: DQRANSPER FUNCTIONS COMPUTET. 
Ceebelbece== PORES, ReSipyuss, AND AZ2n20S (CUMEULED. 
CEmuenies —> CNLY PCLrESeaAND ZEnOS COMBUT ED. 
Cremer = CULT PCLcs AND apsaounsS CONPUTED. 

A 


= 
Sieber CETECN, 1, 2,. 3, 


© 
bu 
PF tas 
® 


eq) 


SGnEeltShGen TRANSFER FUNCTION OPTIONS: 
Grlueoiwie-> NE COMP. TRANSFER FUNCTIONS COMPUTETL. 
melee, —-—- Febbo, RESIDUES, AND ZEROS COMPUTED. 
Cimeroves =~ ChLY POLES AND ZeROS COMPUTED. 
Germor es =-—- CHLY POLES AnD RESLDUES COUPUTED. 


Wiese TA GOURPENSATOR TeAvo SSR ePUNCT ITONSCAN EF 
CCMPUTEL CNLY GF BOTH A REGULATOR 
Sve ena Sti aes Zee 
AND/OR INFUT.) 
cecr. AN CETION;: 1, 2, 3, OR 4. 


J 


Rome FEB D=-FORWARD DISTRIBUTION MATRIX 
wap Rx) «EF INPUT ? 
ieee s HYES'™ OR UNO, 
ILO 
Pomecumrerotnhe LO SYNZTEESEZE A STABLE FELTER (OR REGULATCR) 
EyevecieoerrZeNG THE ORIGINAL SYSTEM? 


(NCTEsWORKS FOR FILTER OR REGULATOR BOUT NOT FCK BCTH 
MISSTHE SALE RUN .} 
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ieee eo Ck TNO. 


t+ 
'o) 


femes Dect Se it ERENT The EVULER-LA 
Peet DECOMPOSITION (@0R CHE 


iereomatoo Ch WYO N, 


° = ~_ a ~~ 
RE S3OG ec 3 ¥\ 4 


Seer ——- COMPUTE TEE PSD OF THE OUTPUTS AND/OR THE 
CC iOLs Ch THE CONTROLLED SYSITEU WHEN FCRCED oY 
PeeersSs END HEASUREMNENI NOTSE. (NOT Ss. Sein kh 
Paeueke OR BND A FILPER MUST Be RESIDENT LN THe 
PaoGhan IC USE THIS OFTION.) 

meee he 2 -—- SANE AS CEFTION 1 AEOVE BUT ONLY PRINT THE 
PeSEUUBS BCE BACH TRANSFER FUNCTION 
USeD een TEE PSD COMPUTATION. 

Meet oOhes -—- NOT LCESIFEDC. 

Seeger AN CETION: 1, 2, OR 3. 


cemen SPSCTSAL LENSITY ~(ESD) OFTION 1 3; 


ea) 


dad 


Bemrece peslke RECULATOR SYNTHESIS ONLY? 
ieee rro CR NON. 


yes 
Pern weep ¢ OF STATES (NS) OF THE SYSTEM MATRIX 
(geet ATRIX). 

Z 
Puree TBE £ OF CCNTROLS (NC) OF THE SYST£M HCDEL 
(1G "= FAT RIX). 

, 


ENTER TEE # OF MEASUREMENTS OR OBSERVATIONS (NO) 
Cont =MATRIL). 


aa) 


Pugeho ihe i OF EROCESS NCISE SOURCES (NG) 


aye) 





feoatee -MATKIX) . 


Peo r aio ER SEITINGS FOR THIS 2UN AnKE AS FCLLCas; 
mueeeme ik {SS I ITF 


welt Ece Pils IP btn 22 LL eEUG 

1 1 C 0 C 1 0 Q 0 0 0 
meee DsTAE FPSD TYU INCEM ITREG NS NC NOB RE 

c Q ‘ C 1 Z 1 @ 0 
MeoER CF SYSTEM = 2 
NOMMEER CF CONTRCLS = 1 
ee 6CF 6CBSERVATICKNS = 2 
BipeerreeereerrOCRSS NCISE SOURCES = QO 


Miveere The SYSTEXR MATRIZ ("F'"-MATRIX) 
EeeMieNomON = # STATES (NS) K # STATES (NS) 
Cpe FIeMENT F( i, 1) 


Mere ELEMENT F( 1, 2)= 


PB ELEMENT F( z, 1)= 


met ELEMENT F( <, 2)= 


Rite SVSTEM MATRIX (SES aMATR ITA) 2. « 


0.0 Herc U0 
=1.CCoO0Cc 0.9 
[LO YCU WISH TO CBANGE THE VALUE OF ANY MATRIX ELEMENT? 
Ries "YES" Ck "NO", 
ho 


Ss) 





Wee N ECCP CY NAM 
er 


Ga 
moet GUL+6 1 Q.0 


Mee or GLGENWALUES. . 3260s. ow UET(SI“F).. 


0.0 , i1-COOOCI+00: 
emreecr RECHT EIGENVECTCEOMATRIN. «2520s aeeee 
U0 = i UNERO Gera Sure 


fecoeeact(D+CO 0.0 


Weerece ee LEST EICENVECTOR MATRIX...... aN ieee 
0.0 1-000CCCD+00 
pie COUCUCCE+CO 0.0 


FNIER IHE MEASUFEMENT SCALING MATRIX ("H"-MATRIX). 
CIMENSION = # OBSERVATIONS (NO) X # STATES (KS) 
THE ELEMENT H( 1, 1)= 


C 
MEO ELEMENT H( 1, 2) = 
0 
IHE ELEMENT H( Z, 1)= 
¢ 
TEE ELEMENT H( Z, 2)= 
2 
-1 
TES MEASUREMENT SCALING MATRIX ("H"-MATRIX)... 
0.0 0.0 
C.0 -1.0CC00 
IC YCU WISH TO CEANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 
no 


apy, 





Perot et UR PM ENT SCALING MATRIX. ..H (BAR) #T.. 
C.0 0 
aeeeeucUgOD+O00 QO. 


Peete OUTPUT MEASUREMENT COST MATRIX (7 A"N-MATRIX). 
PetehorON = # OESERVATIONS (NO) X # OBSERVATICKRS (NC) 
meer here nT A( i, 1)= 


ea) 


c 
TEZ ELEMENT A( 1, 2)= 
0 
fee PEeUnNT OA 2, 1)= 
C 
TEE ELEMENT A( Zz, 2)= 
uy 
THE OWIPUT MEASUREMENT CCST MATRIX ("A"-MATRIX)... 
C.0 0.0 
G0 4.0CC00 
IC YCU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO", 
no 
GEE CT COST MATRIX....-.---- Se ier os: « 
0.0 0.0 
0.0 0.40 00L+01 


Pieeneine, CCNT RCL DISTRIBOEION MATRIX ("G"-MATRIX). 
DMMENSTON = ¢ STATES (NS) X # CONTROLS (NC) 
TEE ELEMENT G( 1, 1) 
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no 


oe ote noe DicthisUuliON MATRIX ("C= aAT 


mee CUSSZSA TO CHANGE TEE VALUE OF ANY MATRIX 


meee ees Ch MNO, 


Meereodes CCONTRCL COST @WEIGHTING MATRIX ("BY-MATRIX) 
PinewotON = + CONTROLS {NC) X # CONTROLS (XC) 


THE ELEMENT B( 1, 1)= 


TEPeONTROL COST MATERIA. «on cw ec Be ae 
1.CC000 


Bomeeecu Wisi TO CEANGE TEE VALUE OF ANY MATRIX ELESENT? 


Ive oats Che"NOi. 
| Meee wTROL DISTEIBUTLION MATRIX. ccc ccc ceeGes 
Omeb 
OeeieGOL+ 0 1 


meme CCNTROL DISTRIBUTION MATRIX.....-11*G.. 
ec Comet D+ 00 
0.0 


THE CENTROL COST MATRIX. ccc ccc we we ccc cee oes 
O. 10COL+01 
CFEN LCCP TRANSFER FUNCTIONS... 


TF FCH IKEUT NOjw 1 AND OUTFUT NO. 1: 


H@ FINGZIE ZEROS. TF GAIN = 0.0 
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Eve Lt zs on So py eS 
S588 (4) IMAG (E) 
Orrcmeyt3( 1.0C000C) { 0.04 EXP (A*T) *COS(B¥*T) 
( 0.C )+3(-1.060000) ( 0.0 ) EXP (A*T) ¥SZN(B¥T) 
TF FCR INEUT NO. 1 BND OUTEUT NO. 2: 
OALER CF X\UMERATOR = 1 TF GAIN = -0.1900D+01 


NUMERATCR EIGENVALUES (INCLULING EXTRANEOUS ZERO VALUES): 
( 0.0 )tJd | 0.0 ) 
( 0.0 eo ( 0.0 ) 


Peeper UES AT THE FOIES: 


Emel: S PEs S Lee 2 OS 
REAL {A) IMAG (B) 
( 0.0 )4+J( 1.000000) ( -1.€00000) EXP (A* I) *COS (E¥T) 
( 0.C )+9(-41.000000) ( 0.0. ) EXP (A¥I) *SIN (E*T) 


feeemene INTHOCE2 ON ELGENVALUE NO; 4 


-1,9€Z36€D+CO 3.464€12D-03 -2.469867D+00 1.508857D+00. 
3.4E€483ED-03 3.762172D-02 -1.491143D+00 2.5001025+00 
-4.4150415-15 -3.208€43D-13 -1.962366D+00 3.621151D-03 
5.2€1648D-11 -1.267812D-17 3.621125D-03 3.762121D-02 


PYG RMNWSYSTEM CF CFEFIIMNAL REGULATOR. «2022 cee 


EULER-LAGRANGE EQUATIONS HAVE A REAL EIGENVALUE 
AT OR NEAF ZERO. 


GomeeermerriNalL REG. E-VALUES...DET(SI-FtG*C) .. 
Ore 2 OY a lero O OCD + OG: 


GmEpGGrehiGHT ZIGENVECTOR MATRIX... cccce sce Meee 
meemueceoiD—C1l 64035 1€5D-01 
poereescvlL-O1 861 .000CCED+00 





PGE a ame Ga 


Peer rCl £liGENVACTCR MATRI 
“ 


=o —~ Be ory ~~ ~ ~~ — 5 ~—; 7 _ e —y == 2 ” — ay 
i tees ie = 5 = ins i oh 1 d ie 
Ee wole oh bes NJ DES be ae Z=y eClLOR Merwe A ea ee ¥ 


C Das. 
MeeebtCO 2.16 7EESDEOO 
meeelru0 2553525E62D+00 


eee GEeliNal PEELTBACK GAIN CCNIROL MATRIX...C=BINV*GI*S... 
= 3 


Sete ool sO CU 0095 70400 


Meme CENT RCTL GAINS... 6. eee ces aces sCELe. 
meomercoc CDtO0 =-5. iC9451D+00 


weet OcrD LCORSLYNAMZCS MATRIX «5.00 «teGeC. . 
0.0 1.000CCCD+00 
See ts WD+ CO —3 .0986°56D+00 


ANALYSIS COMPLETE. CTO YOU WANT ANOTHER KUN? 
TYPE “YES OF "NO", 
yes 
CO YOU WISH TO SAVE THE “E"-MATRIX FROM THE LAST 
BUN TO BE USED IN THE FOLLOWING RUN? 
NCTE: THE MATRIX WILL EE FEDISPLAYED AT 
TEE FRCPER INPUT SEQUENCE INTERVAL 
AND YCU WILL HAVE THE CETION OF CHANGING 
INCTIVILUAL MATRIX ELEMENTS. 
TYPE "YES" OR "NC*, 
yes 
CO YOU WISH TO SAVE THE "“HY-MATRIX FRCM THE LAST 
FUN TO BE USED IN THE FOLLOWING RUN? 
NCTE: THE MATRIY WILL EE FEDISPLAYED AT 
TEE FRCPER INPUT SEQUENCE INTERVAL 
AND *¥CU WILL HAVE THE CETION OF CHANGING 
INCDIVICUAL MATRIX ELEMENTS. 
TYPE "YES" OR "NC", 
yes 
CC YOU WISH TO SAVE THE "G"-MATRIX FROM THE LAST 


Or 
(a) 





q 


cere i moeeD DN He ePCLICWING RUN? 


ere Cee ee ee PEO ES PLAYED AT 
(eqeeeertne THFUL SEQGUECC: LATERY al 

Pee YCU Witi HAW THS CYTION CF CHANGING 
uo twee OA Loe nn oe ENTS. 


ONG! 


be@ run again, 


OL the coSt matrix element 


Since the Same ELrogram options ére 


values, 


“eos 
ed 


With only @ change in ore 


s 
a 
wy ad 


terminal output was deleted up +o tne 


peirnt where program calculations resume 


in crder tc avoid redundancy. 


Onperr OF SYSTEM = z 

NUMETR CF CENTRCLS 1 

MORSE Er CF CESERVATICNS = <z 

NMimiesh Cr EROGESS NCISE SOURCES = Q 


Moe ooo A tha A (TPN = MATROX) << « 
1.00C00 
0.0 


C.0 
-1.CCOCC 


IC YCU WISH TO CEANGE THe VALUE Of ANY MATRIX 


ieee bo eRe "NON, 
EO 
Cee velCeCr OY NANI CS MATRIX< wc csescswcessacealae 
0.0 O.100CI+01 
—OeelGeor +! | 0 .0 


Clin Lecr 


PMEGE NIH INUE Se slele ses cles DO L{ol=r). . 
1.CO00CI+00: 


CEEN 
0.0 


1. 9O00CCD+C0 


eee ctl rr iGeny SCLOR MIATRE AS s. 666s Lou ss 
tee OuCCUDtU0 


0.0 
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PEW ENT 2 








0.0 
ie VOMmGtL+CQ G 


no 


Gee eee ee TCE VECTOR MATRixX«. ees Te iNy.. 


1900 0CCD+00 
20 


Peso SCALING HALT RIA ("HY=K4ATRIX) «.. 

GE 0 0.0 

Ga Sree eu UY 

Pee eo Visoe tO CRANGE Eder YVRLUE OF ANY MATPITX ELEMENT? 
fees yes CR NON. 


eerste as «6CSCAL TING PYWAT RIX. «ccc aceancces Gee ee 
0.0 
-0.100CI+01 


NCTAL MEASUREMENT SCALING MATRIX.~.-.H (BAR) *T.. 
¢.0 0.¢ 
Siem goto O. C 


Soe Cemetery) URASU EERE NT COST MATRIA ("A™ “0 ATRIX) « 
ieee > £4 OESERVATIONS (NO) X # OBSERVATICNS (30) 
TEE ELEMENT A( 1, 1452 


C 

Bee Hire NT A( 1, 2)= 
0 

Ret EEEIENT A( 2, 1)= 
C 

Meow NT A(z, 2)= 
4.000 1 


Meee UrrOr sUSASUBBUENT CCST MATRIX ("AN-MATRIX)... 

C.0 0.0 

C.0 4.0C010 

ECmpoGlmenGStelO CFANGE TEE VALUE Of ANY MATRIX ELEMENT? 


Ce 





Gay tes CR NOM, 


rec 
CUE EEE GOSo 1S Mano ot ae eo 
oe 0.0 
0.0 0.40 0CL+01 
TEE CCNTROL DISTRIBUTION MATRIX ("GU-MATRIX)... 
0.0 
1.CCOOC 
[C YCU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NON, 
ho 
ENTER IHE CCNTRCI CCST WEIGHTING MATRIX ("S"-MATRIX) 
CIMENSION = # CONTROLS (NC) X # CONTROLS (XC) 
THE ELEMENT B( 1, 1)= 
; 
MEBCOMMROL COST MANNIX. ....... +B. se 
1.€6000 
IC YCU WISH TO CHANGE THE VALUS OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO". | 
NO 


eee wCCnue nol DISTEIBUTION UATAIX. . «22102 ee Gee 
0.6 


Cee oeureu"| 


meme COnmnolL DESTalsuUrTiIon WATRIX......-TiL*G.. 
limeeooe Cp+00 
0.0 


ereMEERO IL 6GOo | LATRI A... 6 cc occ co ew cc ew ee ole os 
QO. 1CCOL+01 


Creve LOCr TRANSFER FUNCTIONS. .. 
Meee nul NO. ean OUREUT NO. 1: 
NG@wEONeUeegrnCs. Ir GAIN = 0.0 
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re @eaaG (5) 
meee jeu ( 1.0000CC) (ee. 0.) Een CC oes) 
2:6 Vee oe oc C) (oe 0 (} BAe (ae jes 8 (es T ) 


Shien Cr SUAERATCR = 1 TF GAIN = -90.10000+01 


memmereren EIGENVALUES (INCLULING EXTRANEOUS ZGcRO VALUES}: 
( QO.C )+Jd { 0.0 ) 


Meee UES Pe Vit POLES: 


6 ie es REESE ES 
REAL (a) IMAG (E) 
Weeeeo. )+5t 1.CG00TC) ( —-1.000000) EXE (A*T) *COS {B*T) 
( 0.C )+#J3¢-1.0000CC) (080.0 ) EXP (A*T) *SIN (B4T) 
FIGEXSYSTEM CF CETIMAL FEGULATOR......-. ae 


mec ce 1 h+0C2—9.950 1206-0 1; 


GC=1CCPE FIGHT £IGENVECTOR MATHIA.. 2-2 siLisle «| 6 
feo ese ~ Ot —~7/.088 723D-01 
moeececrzep=O01 7 .053368D-01 
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In itS present interactive form, OPTSYSX has been trans- 
ferred fren the arena of high-level applied mathematics and 
Rumericai analysis tc the level cf control system angineer 
and students. It now represents an ¢€¢ven more rowerful 
educational tool, étle to rapidly and effectively unicck 
Tany agisunderstcod linear systems mathematical relaticn- 
ships. 

em wamerteetiee SValuacecmeetethe computational arilities 
SemOrIsiSX, the procram was tested using an 82 X 8&2 matrix 
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experimental Fichter aircraft derivative, provided by 
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Frogram results exhipited perfect eigenvalue correlation 
with these cbhtained from the John Edwards Centrol Freqran. 
Additicnally, OPTSYSX provided complate longitudinal mecdai 
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Gueeeeecne NCn-IKSL library of subroutines, makirg it 
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